Aperiodicity-Induced Second-Order Phase Transition 
in the 8-State Potts Model 
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We investigate the critical behavior of the two-dimensional 8-state Potts model with an aperiodic 
distribution of the exchange interactions between nearest-neighbor rows. The model is studied 
numerically through intensive Monte Carlo simulations using the Swendsen-Wang cluster algorithm. 
The transition point is located through duality relations, and the critical behavior is investigated 
using FSS techniques at criticality. For strong enough fluctuations of the aperiodic sequence under 
consideration, a second order phase transition is found. The exponents jS/u and 7/1^ are obtained 
at the new fixed point. 



The study of the influence of bond randomness on 
phase transitions is a quite active field of research, mo- 
tivated by the importance of disorder in real experi- 
ments 0] . According to the Harris criterion , quenched 
randomness is a relevant perturbation at a second order 
critical point when the specific heat exponent a of the 
pure system is positive. 

The analogous situation when the pure system exhibits 
a first order phase transition was also studied. Imry and 
Wortis first argued that quenched disorder could induce 
a second order phase transition |^] , and it was shown that 
in two dimensions, an infinitesimal amount of randomly 
distributed quenched impurities changes the transition 
into a second order one |4j. The first large-scale Monte 
Carlo study of the effect of disorder at a temperature- 
driven first order phase transition is due to Chen, Fer- 
renberg, and Landau. These authors studied the 2D 8- 
state Potts model (which is known to exhibit a first order 
phase transition when the number of states q is larger 
than 4 |^). They first showed that the transition be- 
comes second order in the presence of bond randomness, 
and obtained critical exponents very close to those of the 
pure 2D Ising model at the new critical point |6| . On the 
other hand, drastically different results were obtained for 
random lattices 

The essential properties of random systems are gov- 
erned by disorder fiuctuations. All physical quantities 
depend on the configuration of disorder, and the study 
of the influence of randomness requires an average over 
disorder realisations. Among the systems where the pres- 
ence of fiuctuations is also of primary importance, ape- 
riodic systems have been of considerable interest since 
the discovery of quasicrystals ||] . They are built in a de- 
terministic way, making any configurational average use- 
less, and their critical properties have been intensively 
studied (for a review, see Ref. ||]). In layered systems, 
aperiodic distributions of the exchange interactions be- 
tween successive layers in the Ising model have been con- 
sidered ||l^, leading to unchanged universal behavior or 



to modified critical properties, depending on the aperi- 
odic series under consideration. The major result was 
obtained when Luck, generalizing the Harris criterion to 
layered perturbations, proposed a relevance criterion for 
the fiuctuating interactions . According to Luck's cri- 
terion, aperiodic modulations may be relevant, marginal, 
or irrelevant, depending on the correlation length expo- 
nent v of the unperturbed system and on a wandering 
exponent co which characterizes the fiuctuations of the 
couplings around their average |l2|. Systematic studies 
of the critical properties for irrelevant, marginal, and rel- 
evant aperiodic perturbations have then been achieved in 
the extreme anisotropic limit p3| . 

In this letter, we report results of Monte Carlo simu- 
lations of the two-dimensional 8-state Potts model with 
an aperiodic modulation of exchange couplings between 
nearest-neighbor layers. Our aim is to study the effect of 
such a distribution on the nature of the phase transition. 
In particular, we ask if the fiuctuations are able to induce 
a second order phase transition. The Hamiltonian of the 
system with aperiodic interactions can be written 



(1) 



where the spins cr^, located at sites i, can take the values 
a = 1, 2, . . . , g and the sum goes over nearest- neighbor 
pairs. The coupling strengths are allowed to take two 
different values Kq = K and Ki — Kr. They are dis- 
tributed according to a layered structure i.e. the dis- 
tribution is translation invariant in one lattice direction, 
and follows an aperiodic modulation {/fc} of digits fk — 
or 1 in the other direction: In layer k, both horizontal 
and vertical couplings take the same value Kr^'' (Fig. 
The sequence of digits {fk} is generated through iteration 
of substitution rules. The Thue- Morse (TM) sequence is 
obtained by substitutions on digits: '5(0) = 01, 
1 — * S{1) = 10, while the so-called paper folding se- 
quence (PF) is generated through substitutions on pairs 
of digits: 11 ^ 5(11) = 1101, 10 ^ 5(10) = 1100, 
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FIG. 1. Layered aperiodic modulation of the coupling 
strengths on the square lattice and dual system. 

Our particular choice of coupling distribution makes 
it possible to determine exactly the critical point by du- 
ality arguments. Consider a system of L layers with a 
distribution {/fc}, made from a succession of vertical- 
horizontal (V-H) bonds when read from left to right 
(Fig. and let us write its singular free energy den- 
sity fs{KQ,Ki]{fk})- Under a duality transformation, 
the strong and weak couplings Ki are replaced by weak 
and strong dual couplings Ki respectively. Since a verti- 
cal bond on the original lattice becomes horizontal on 
the dual system, the same V-H bond configuration is 
recovered for the transformed system when the distri- 
bution is read from right to left. One thus gets the 
same type of system, but a reverse distribution 
so that the free energies of the two systems are the 

same: /,(Xo, ^i; {A}) = /.(i^o, i^i; {/L+i-fc})- The 
sequences considered here have the interesting prop- 
erty that the reverse distribution corresponds to the 



original one if one interchanges perturbed and unper- 
turbed couplings Ki ^ Kq: fs{Ko,Ki]{fL+i-k}) = 
fs{Ki,Ko',{fk})- The system being self-dual the criti- 
cal point, if unique, is exactly given by the critical line 
Kqc — Kic of the usual anisotropic model |14]: 
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l)(e^^'--l)=g. 



(4) 



One should mention there that the required symmetry 
property of the sequences holds for odd iterations in the 
case of TM, and works in the case of PF if one omits the 
last digit, which simply introduces an irrelevant surface 
effect in the simulations. 

We performed extensive simulations oi L x 2L lat- 
tices (16 < 2L < 512) with periodic boundary condi- 
tions in one direction (vertical) and free boundaries in the 
other {2L columns). The Swendsen-Wang cluster flipping 
method jl^ was used. Between 2 x 10^ (smaller lattice 
sizes) and 6 x 10^ (larger lattice sizes) Monte Carlo steps 
(MCS) per spin were performed (this is always larger 
than 10^ times the correlation time, and seems sufficient 
in order to produce reliable thermal averages). 

The order parameter is defined by the majority orien- 
tation of the spins n& : 



M = (m) = 
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(5) 



Here, pmax = {maxa-ipa-)), where pa is the density of 
spins in the state a and (. . .) denotes the thermal aver- 
age over the Monte Carlo iterations. The susceptibility 
is given by x = KV{{m?) — (m)^). Although local or- 
dering mechanisms are not yet clarified in aperiodic sys- 
tems, we expect a unique transition temperature for all 
the columns, so we used average quantities in order to 
reduce fluctuations. 

The flrst task is to identify the order of the transi- 
tion. For that purpose, we made some preliminary runs 
at several temperatures first of all in order to confirm nu- 
merically the location of the critical point, and then to 
have a general picture of the phase transition. The ex- 
amination of the energy autocorrelation time shows that 
it is diverging in the case of TM sequence, while it seems 
to remain more or less bounded for PF. It is consistent 
with a first order transition in the first case, and a sec- 
ond order one for the latter. We have further estimated 
temperature-dependent effective exponents for the aver- 
age magnetization and suceptibility. This can be done by 
comparing the data at two different sizes L and L' = L/2: 
Assuming the following scaling form for the average mag- 
netization ML{t) = L-P/^MiLr), where t - K^l 
we define the quantity p7[ 



\n Ml/Ml 
In L/L' 



(6) 



Close to Kc, this can be expanded in powers of Lf^ , lead- 
ing to 
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which defines an effective exponent. As the critical point 
is approached and in the thermodynamic hmit, it evolves 
towards — ^. The analogous quantity can be computed 
for the susceptibility. The results are shown for TM on 
Fig. i 
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FIG. 2. Temperature-dependent effective exponents for 
the average magnetization and suceptibility (TM sequence) 

In 



estimated from the ratios 
YL(t) = '""^/f^- (7^1/7^0 



XL[t) = -"-r"^' and 
- r = 5). The sizes vary from 
21/ = 16 to 2L = 256, and the arrows indicate the crossings 
between successive sizes which give rough estimates for the 
critical exponents. 

The successive estimates oi ^ = d — yt and ^ — 2yh — d 
clearly evolve towards the values and 2, characteristic 
of a first order phase transition. The scaling dimensions 
associated to the temperature and magnetic field, yt and 
yh, indeed take a special value equal to the dimension d 
of the system In the case of the PF sequence, the 
behavior is drastically different, and this first analysis 
does not allow any conclusion. 

Once the qualitative description of the phase transi- 
tion was made, our strategy was to use finite-size scaling 
(FSS) techniques in order to get more accurate results. 
We made runs for systems of larger sizes, and in a L x 4L 
geometry (4L in the aperiodic direction, going from 8 to 
1024), for which we estimated the number of MCS/spin 
from the preliminary runs. We have moreover studied a 
periodic system (PS) with alternate couplings Kq and Ki 
(i.e. the same critical point given in Eq. ^), and which 
is a "first order reference" system. This is illustrated in 
Fig. 1^ where the energy autocorrelation time r is plot- 
ted for the three samples. In the case of TM, r diverges 
exponentially as expected for a first order phase transi- 
tion, although it is always quite small compared to the 
periodic system. For PF, the data are compatible with 
a power law with a very small dynamical exponent as 



expected for cluster algorithm simulations at a second 
order phase transition. 
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FIG. 3. Energy autocorrelation time r at K'^ {r — 5). For 
TM, the dashed line is a fit to an exponential behavior. 
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FIG. 4. Log-log plots of M and x vs L for a periodic ref- 
erence, and for TM and PF sequences {r — 5). For M, error 
bars are smaller than the symbol sizes. 

The crude data for M and x furthermore show that, 
in the case of the TM sequence, a cross-over appears 
between small sizes where the data more or less follow 
the same behavior than PF, and large sizes where the 
first order regime analogous to PS is well established 
(Fig. ^). Certainly, a careful procedure is needed for 
a reliable determination of the critical exponents. From 
the log-log curves between ALmin and '^Lmax = 1024, one 
determines an effective exponent x{Lmin) for each quan- 
tity, then the smaller size is cancelled from the data and 
the whole procedure is repeated until the three largest 
sizes only remain. The effective exponent is then plot- 
ted against (Fig. |^). The critical exponent is finally 
deduced from the extrapolation at infinite size. The nu- 
merical results are given in Table ^. All of them are in 
agreement with the scaling law d = 2(3 /v + ^ /v, within 
the precision of the results. 
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To summarize, we have shown from numerical simu- 
lations that the fluctuations introduced by an aperiodic 
modulation of exchange interactions is liable to induce 
a second order phase transition in a system which origi- 
nally exhibits a first order transition. 
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FIG. 5. Effective size-dependent exponents for the two ape- 
riodic sequences (r = 5). 

TABLE I. Exponents associated to the magnetization and 
susceptibility for the three examples considered in the text. 
The uncertainties given there are rough estimations corre- 
sponding to the standard deviation for the fit of the data in 
the whole range of sizes. 





PS 


TM 


PF 


13/u 


0.05 ±0.03 


0.05 ±0.04 


0.48 ±0.03 




1.99 ±0.08 


1.96 ± 0.06 


1.01 ±0.04 



From Monte Carlo simulations, we have strong evi- 
dences in favor of a second order regime for PF sequence. 
This type of effect was already known since the work of 
Imry and Wortis in the case of a random distribution 
where the fluctuations are unbounded and can be charac- 
terized by a wandering exponent cji-and = 1/2. Here, the 
same type of behavior is induced by a smoother perturba- 
tion, namely by the PF sequence which exhibits only log- 
arithmic fluctuations (to = 0) while the bounded fluctu- 
ations generated by the TM sequence {uj — -co) are not 
strong enough to destroy the first order transition. The 
same type of problem on a quasicrystal is currently un- 
der investigation by Ledue et al and the transition seems, 
in this case also, to remain of first order We may 

thus infer that Luck's criterion can probably be applied 
to first order phase transition. Here, we can replace 
by 1/cZ p^ , ^0| in the criterion in order to compare ther- 
mal fluctuations to those introduced by the distribution 
of couplings. Luck's cross-over exponent then becomes 
(f> = 1 + (lu ~ 1) / d, and the aperiodicity can induce a 
second order phase transition when cj) > 0. This is in 
agreement with the results of our simulations. We can 
finally mention that the study of local order parameter 



is currently under investigation. 
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